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Abstract 


We  calculate  the  change  in  the  density  of  states  due  to  a  biased 
resonant  tunneling  structure.  The  maximum  of  the  density  of  states  near 
resonance  gets  shifted  towards  low-energy  side  compared  to  the  unbiased  case, 
as  is  the  transmission  coefficient,  although  the  two  need  not  be  identical. 
For  the  case  of  asymmetric  barrier  heights,  the  left- right  symmetry  of  the 
density  of  states  is  broken  when  the  field  is  non- vanishing . 
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In  recenc  years  sciencists  nave  become  increasingly  aware  of  the 
importance  of  resonant  tunneling  structures  (RTSs)  in  electronic  and  opto¬ 
electronic  device  applications  and  their  possible  role  m  fundamental 
advances.^  ^  Therefore,  many  studies  of  static  and  dynamic  aspects  of 
resonant  tunneling  structures  have  been  undertaken.  A  common  static  quantity 
of  interest  is  the  transmission  coefficient  T(E);  a  related  physical  quantity 
is  the  width  of  the  resonance  peak,  which  is  inversely  proportional  to  the 
lifetime  of  the  resonant  state.  A  second  important  static  quantity  is  the 
density  of  states  ,  a  knowledge  of  which  is  essential  in  understanding 
transition  probabilities,  dielectric  functions  and  absorption  and  luminescence 
characteristics . 

Recently  the  local  density  of  states. 

N(E..x)  -  ^  5(E-E^)  |0^(x)  1^  ,  (1) 

n 

in  a  RTS  has  been  obtained  and  analyzed  in  various  limits.^ Here  E  are  the 

energy  eigenvalues  of  the  system  and  are  the  corresponding  eigenstates. 

The  global  density  of  states  obtained  by  integrating  Eq .  (1)  following  Ref.  5 

would  be  identically  zero.  We  therefore  follow  a  different  scheme  to 

a 

calculate  the  global  density  of  states, 

N(E)  -  ^  5(E-E^)  (2) 

n 

In  a  big  box  of  size  L  (L  -*  «>)  ,  the  energy  levels  form  a  quas  i  -  continuous 
spectrum.  Introduction  of  a  structure  inside  the  box  changes  the  spacing 
bc^’-foen  Che  levels  and  produces  a  change  AN(E)  in  Che  global  density  of 


states,  which  we  calculate.  In  the  neighborhood  of  a  resonance,  the  changes 
in  the  spacing  of  the  energy  levels  produces  a  pronounced  change  in  density  of 
states.  If  the  resonances  are  sufficiently  narrow,  AN(E)  and  T(E)  coincide; 
however,  for  broad  resonances  this  is  not  true. 

For  the  sake  of  completeness,  va  briefly  review  our  method  for  obtaining 

g 

the  density  of  states  AN(E)  for  an  unbiased  RTS.  For  a  flat  box  extending 

from  X  -  0  to  X  -  L,  the  density  of  states  in  k-space  is  N^(k)  -  L/jt,  where  E 

2  2 

-  h  K  /(2m*)  is  the  energy  and  m*  the  effective  mass.  Suppose  now  that  the 
RTS  is  placed  in  the  middle  of  the  box,  at  x  -  x.j^  »  L/2 ,  thereby  squeezing 
more  states  into  some  energy  region  and  depleting  states  in  some  ocher  region. 
Suppose  chat  the  energy  eigenvalues  of  the  system  are  obtained  from  the 
condition  D(k)  -  0  where  D(k)  is  a  determinental  function  of  the  solutions  to 
the  Schrbdinger  equation  (see  Eq.  (9)  below).  The  change  in  density  of  states 
associated  with  the  n-ch  level  having  spacing  Ak^  is 

'•  n  ' 

The  spacing  Ak^  is  easily  obtained  by  finding  the  roots  of  the  eigenvalue 
condition  D(k)  -  0  with  the  use  of  a  Newton- Raphson  method  or  any  ocher 

g 

appropriate  scheme.  As  the  previous  work  in  our  group  emphasized,  the  shifts 
of  the  energy  levels  depend  sensitively  on  the  phase  of  the  wave  function  at 
Che  position  where  the  structure  is  introduced  (i.e.,  at  x^  a  L/2  in  this 
case).  This  phase  dependence  produces  apparently  irregular  spacings  of  the 
levels,  and  one  has  to  calculate  two  "sub-densities"  (since  the  RTS  is  in  the 
middle)  in  the  manner  indicated  above  and  add  the  two  to  obtain  the  total 
density.  As  expected,  for  a  biased  RTS,  calculation  of  AN(k)  and  hence  AN(E) 


is  more  compLicated,  and  Che  phase  of  the  wave  function  in  the  neighborhood  of 


the  structure  has  to  be  constantly  adjusted  (by  changing  the  position  to 
x((E))  CO  get  the  sub-densities  correctly.  The  position  x^  at  which  the 
biased  RTS  should  be  placed  is  obtained  from  the  relation 

k'x^  -  ,  (4) 

where  k  -  (2m*E/]<'^)  ^ ,  x-  -  L/2  and  k'  -  [2m*(E+V  Here  V  is  the 

i  0  0 

potential  drop  across  Che  double  barrier  structure  and  is  taken  to  be  10  meV 
throughout  this  work. 

In  Fig.  1  we  show  the  geometry  of  the  device.  The  RTS  has  an  extension 
-  a^  +  a^  +  d,  where  a^  and  a„  are  the  barrier  widths  and  d  is  the  well 
width.  The  barrier  heights  are  taken  to  be  and  V^.  The  electric  field  in 
the  structure  is  uniform  and  is  F  -  .  The  electric  potential  at  any 

point  X  inside  the  structure  is 

V^(x)  -  V^(x^-x)/x^  (5) 

We  label  Che  regions  of  piecewise  continuous  potential  profiles  by  integers, 
from  0  to  4,  as  indicated  in  the  figure,  and  solve  the  s tationary- state 
Schrbdinger  equation  for  the  envelope  function  in  the  effective  mass 
approximation  where  we  assume,  for  simplicity,  Che  same  m*  -  0.067  m^  (m^  - 
electron  mass)  throughout  the  structure.  In  region  1,  for  example,  the 
Schrbdinger  equation  is 

dVdx^  -  (2m*/M^)  [Vp(x)  +  (V^-E)l  vi(x)  -  0  .  (6) 


This  can  be  reduced  to 


d^i/i(p)/dp^  -  (>v(p)  -  0 


(7) 


whose  solutions  are  the  Airy  fiinctions,  Aj^(/j(x)],  and  the  complementary  Airy 


functions  B^[/3(x)], 


i>(p)  -  A^A^(p)  + 


(8) 


wn 


ere  A^  and  are  two  arbitrary  constants,  p(x)  -  ai(x, -x)  +  (V  -£)x-/V  ] 
1.  L  '  1.  L  3  O  ' 

2  1/3 

and  a  -  [2m*V^/()(  x^)  ]  The  solutions  in  all  the  five  regions  can  be 

obtained  in  a  similar  fashion.  The  eigenvalue  condition  is  the  condition  of 
the  vanishing  of  the  wave  function  at  x  -  L,  so  that 


D(k)  -  A^sink' (L-x^)  +  B^cosk' (L-x^)  -  0 


(9) 


where  x_  -  x,  x- ,  and  A.  and  B,  are  obtained  by  demandinz  the  usual 

2  1  3  u  4.  ■'  ^ 

continuity/  of  the  wave  function  and  its  first  derivative  with  respect  to  x 
across  the  interfaces: 


[A^  ^4!^  ”  ^  [sinkx^j^  kcoskx.j^]^ 


The  2x2  matrix  M  in  the  equation  above  is 


(10) 


(11) 


6 


The  subscripts  on  the  matrices  indicate  the  two  regions  they  connect,  and  L 
and  R  stand,  respectively,  for  the  left  and  right  sides  of  the  interface. 

Equations  (3),  (4),  (10)  and  (11)  enable  us  to  obtain  AN(E)  for  a  biased 
RTS,  which  can  then  be  compared  with  the  transmission  coefficient  T(E) 
obtained  in  the  usual  way.  We  call  the  structure  shown  in  Fig.  1  as  a  tilted 
box  with  a  structure  (TBWS).  We  define  a  background  potential  profile  called 
the  tilted  box  (TB)  for  which  -  0  in  Fig.  1.  The  difference  in 
density  of  states  between  TBWS  and  TB  gives  AN(E),  which  can  be  compared 
directly  with  T(E). 

Figures  2-4  show  T(E)  and  AN(E)  for  a  double  barrier  device  of  barrier 
widths  50  A  each  and  heights  200  meV  each,  and  a  well  width  of  100  A.  Figures 
2  and  3  show  T(E)  and  AN(E)  for  the  first  two  bound  state  resonances,  whereas 
Fig.  4  is  for  energies  above  the  barrier  energy.  A  comparison  with  the 
unbiased  case  shows  chat  both  T(E)  and  AN(E)  get  shifted  to  lower  energies 
chan  the  corresponding  unbiased  case. 

For  the  sake  of  completeness,  we  comment  on  the  physical  origin  of  the 
energy  shift  of  transmission  resonances  due  to  the  electric  field.  The 
electron  resonance  energy  ^s  a  compromise  between  Che  increased  kinetic  energy 
due  to  a  spatially-varying  potential  and  the  lowering  of  the  potential  energy 
brought  about  by  the  field.  The  electron  wave  function  wiggles  a  lot  to 
accommodate  the  increased  kinetic  energy  and  lowers  its  total  energy.  The 
electron  resonance  energy  decreases  linearly  for  both  the  ground  state  and  the 
first  excited  state  for  the  electric  field  considered  by  us.  The  magnitude  of 
the  rate  of  decrease  with  the  field  is  larger  for  the  ground  state  than  the 


first  excited  state. 


Finally,  in  Fig.  5  we  show  our  results  for  an  asymmetric  double  barrier 
structure  without  a  field  (V^  =  0)  and  with  a  field  (V^  =  10  raeV)  .  The  dot- 
dashed  curves  in  Fig.  (7a)  and  (7c)  represent  the  field-free  case,  and  these 


have  left-right  symmetry  with  respect  to  the  interchange  of  the  two  barriers. 
In  the  presence  of  a  field,  this  symmetry  is  broken.  The  dashed  lines 
correspond  to  the  case  -  100  meV  and  -  200  meV,  whereas  the  dotted  line 
corresponds  to  the  permuted  case  -  200  meV  and  -  100  meV .  This 
asymmetry  is  important  in  calculating  tunneling  currents  and  transition 
probabilities . 
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Figure  Captions 

1.  inetry  of  the  tilted  box  with  a  structure  (TBVS)  .  The  box  extends  froni 
X  -  0  to  X  -  L  (we  take  L  to  be  of  the  order  of  1.5  x  10°  A  to  1.0  x  10^  A 

in  this  work)  whereas  the  structure  extends  from  x  =  x.,tox--x  rx, 

1  13 

where  x.  -  a.,  +  a.,  +  d.  Tie  external  electric  field  is  F  -  V  /x  We 

-  o  3 

cake  -  10  meV  throughout  this  paper.  The  zero  of  the  energy  for  this 
and  the  following  figures  (Figs.  2-5)  is  taken  to  be  the  bottom  of  the 
left-most  part  of  the  box. 

2.  Transmission  coefficient  T(E)  and  density  of  states  AN(E)  for  a  symmetric 

double  barrier  structure  (DBS)  in  an  anplied  electric  field.  The  barriers 

are  each  50  A  wide  and  200  meV  high,  and  the  well  is  100  A  wide.  The 

energy  range  shown  is  in  the  neighborhood  of  the  first  resonance  energy. 

The  middle  panel  shows  AN(E)  for  a  tilted  box  (TB) ,  (solid  curve)  and  a 

TBWS  (dashed  curve). 

3.  T(E)  and  AN(E)  as  in  Fig.  2,  but  for  the  second  resonant  state. 

4.  T(E)  and  AN(E)  as  in  Fig.  2,  but  for  the  energies  above  the  barrier 
energy . 

5.  Transmission  coefficient  T(E)  and  density  of  states  AN  (  E )  for  an 
asymmetric  DBS.  The  barriers  and  the  well  are  each  50  A  wide,  and  the 
barrier  heights  are  100  meV  and  200  meV.  The  dash-dotted  curves  show  T(E) 


and  AN(E)  without  the  field,  and  these  curves  exhibit  left  and  right 
degeneracies.  The  dotted  curves  are  for  -  200  raeV  and  -  100  meV , 
whereas  the  dashed  curves  are  for  -  100  meV  and  -  200  meV.  The 
solid  curve  in  the  middle  panel  is  for  a  tilted  box,  as  explained  in  the 
caption  to  Fig.  2. 
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